Quantum phase transitions from Solids to Supersolids in bi-partite lattices 
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We studied some phases and phase transitions in an extended boson Hubbard model slightly 
away from half filling on bipartite lattices such as honeycomb and square lattice. We find that 
in the insulating side, different kinds of supersolids are generic stable states slightly away from 
half filling. We propose a new kind of supersolid: valence bond supersolid. We show that the 
quantum phase transitions from solids to supersolids driven by a chemical potential are in the same 
universality class as that from a Mott insulator to a superfluid, therefore have exact exponents 
z = 2, v = 1/2, 77 = ( with logarithmic corrections ). Comparisons with previous quantum Monte- 
Carlo (QMC) simulations on some microscopic models on a square lattice are made. Implications 
on possible future QMC simulations are given. 



1. Introduction: A supersolid is a state with both su- 
perfluid and solid order. Recently, by using the torsional 
oscillator measurement, a PSU's group lead by Chan ob- 
served a marked 1 ~ 2% Non-Classical Rotational Iner- 
tial ( NCRI) even in bulk solid 4 He at - 0.2K [1]. The 
NCRI is a low temperature reduction in the rotational 
moment of inertia due to the superfluid component of the 
state [2j. If this experimental observation indicates the 
existence of 4 He supersolid remains controversial [H, 0. 
However, it was established by spin wave expansion [5[ 
and quantum Monte-carlo (QMC) @, 0, HI simulations 
that a supersolid state could exist in an extended boson 
Hubbard model (EBHM) with suitable lattice structures, 
filling factors, interaction ranges and strengths. But so 
far, the universality classes of the quantum phase transi- 
tions from solids to supersolids have never been studied. 
In this letter, by using the dual vortex method devel- 
oped in we investigate some phases, especially su- 
persolids and quantum phase transitions in an extended 
boson Hubbard model (EBHM) on bipartite lattices such 
as honeycome and square lattice near half filling. Al- 
though the SS in lattice models is different from that in 
a continuous systems, the results achieved in this paper 
on lattice supersolids may still shed some lights on the 
possible microscopic mechanism and phenomenological 
Ginsburg-Landau theory of the possible He supersolids 
0| . The EBHM in honeycomb and square lattices could 
be realized in ultracold atoms loaded on optical lattices. 
So the results achieved in this paper may have direct im- 
pacts on the atomic experiments. 

The EBHM with various kinds of interactions, on all 
kinds of lattices and at different filling factors is described 
by the following Hamiltonian [Io| : 

H = -tJ2 ( b i b j + h - c -) - n * + J n ^ Ul - 



Vi UiUj + V 2 ^2 UiUk 

<<ik» 



(1) 



where n, = b\bi is the boson density, t is the nearest 
neighbor hopping amplitude. U,Vi,Vz are onsite, near- 



est neighbor (nn) and next nearest neighbor (nnn) inter- 
actions respectively, the • • • may include further neigh- 
bor interactions and possible ring-exchange interactions. 
A supersolid is defined as the simultaneous orderings of 
ferromagnet in the XY component and CDW in the Z 
component. Honeycomb lattice is not a Bravais lattice, 
so may show some different properties than those in a 
square lattice. 

2. The dual vortex method. The EqnQ] with only 
the onsite interaction was first studied in Ref . [Io| . The 
effects of long range Coulomb interactions on the transi- 
tion was studied in ■ Very recently, the most general 
cases EqnQ] in a square lattice at generic commensurate 
filling factors / = p/q (p,q are relative prime numbers ) 
were systematically studied in [l3l ] . After performing the 
charge- vortex duality transformation, the authors in (l3j 
obtained a dual theory of EqnQ]in term of the interacting 
vortices ipi hopping on the dual lattice subject to a fluctu- 
ating dual " magnetic field" . The average strength of the 
dual " magnetic field " through a dual plaquette is equal 
to the boson density f = p/q. This is similar to the Hof- 
stadter problem of electrons moving in a crystal lattice 
in the presence of a magnetic field. The magnetic space 
group (MSG) in the presence of this dual magnetic field 
dictates that there are at least g-fold degenerate minima 
in the mean field energy spectrum. The q minima can 
be labeled as tpi, I = 0, 1, • • • , q — 1 which forms a q di- 
mensional representation of the MSG. In the continuum 
limit, the final effective theory describing the supercon- 
ductor to the insulator transition in terms of these q order 
parameters should be invariant under this MSG. In this 
letter, we will extend the dual vortex method to study 
the EBHM EqnQ] in honeycomb lattice at and slightly 
away from q = 2 (Fig. la) 

The dual approach is a MSG symmetry-based ap- 
proach which can be used to classify some phases and 
phase transitions. But the question if a particular phase 
will appear or not as a ground state can not be addressed 
in this approach, because it depends on the specific values 
of all the possible parameters in the EBHM in EqnQ] So 
a microscopic approach such as Quantum Monte-Carlo 
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(QMC) may be needed to compare with the dual field 
theoretical approach. The dual approach can guide the 
QMC to search for particular phases and phase transi- 
tions in a specific model. Finite size scalings in QMC 
can be used to confirm the universality class discovered 
by the dual approach. 

3. The effective action and order parameters in the 
dual vortex picture. The dual lattice of the honeycomb 
lattice is a triangular lattice. Two basis vectors of a prim- 
itive unit cell of the triangular lattice can be chosen as 
ai = x, a% = — ix+ -Tpjy, 3d = ai + 3 2 as shown in Fig. la. 
The reciprocal lattice of a triangular is also a triangular 
lattice and spanned by two basis vectors k = fci&i + k 2 b 2 
with 61 = x + -^=, b 2 = -^y satisfying b t ■ 3j = The 
point group of a triangular lattice is Cq v ~ Dq which 
contains 12 elements. The two generators can be chosen 
as Ce = R n /s,Ii. The space group also includes the two 
translation operators along 3\ and 3 2 directions T\ and 
T2. The 3 translation operators Ti,T2, T^, the rotation 
operator R n /3, the 3 reflection operators Ii,I 2 ,I ( i, the 
two rotation operators around the direct lattice points A 
and B: i^/g, #11/3 °^ ^ e MSG are worked out in [lit - 
It can be shown that they all commute with TL V . How- 
ever, they do not commute with each other, for example, 
T{T 2 = toT d , T X T 2 = w 2 T 2 Ti where w = e i2w/ . 

In the following, we focus on q — 2 case [l5[ where there 
is only one dual vortex band E(k) — — 2i(cos fci+cos k 2 — 
cos(fci + fc 2 )). Obviously, E(k\,k 2 ) = E{— ki, — fc 2 ) = 
E{k 2l k\). There are two minima at k± — ±(7r/3, 7r/3). 
Let's label the two eigenmodes at the two minima as 
ipa/b- How the two fields transform under the MSG was 
derived in [l5j |. 

Moving slightly away from half filling / = 1/2 cor- 
responds to adding a small mean dual magnetic field 
H ~ Sf = f — 1/2 in the action. It can be shown that 
inside the SF phase , the most general action invariant 
under all the MSG transformations upto quartic terms is 



C S f = K fl f . - iA^cf +r\ip a \ 2 + -L^a^Ax 

a— a/b 

+ 7o(iv^i 2 + iv^i 2 ) 2 - 7i(^ai 2 - m 2 ) 2 + ■■■ 



where A^ is a non-compact U (1) gauge field, e is a dimen- 
sionless coupling constant depending on t, U, Vi, V 2 ■ ■ ■ in 
EqnJT] Upto the quartic level, with correspondingly de- 
fined ip a /b in a square lattice, Eqn[2] is the same as that 
in the square lattice derived in [131 ] . 

Because the duality transformation is a non-local 
transformation, the relations between the phenomeno- 
logical parameters in Eqn|2] and the microscopic param- 
eters in Eqn[T] are highly non-local and not known. For- 
tunately, we are still able to classify some phases and 
phase transitions and make some very sharp predictions 
from EqnJ5] without knowing these relations. If r > 0, 



the system is in the superfluid state < -0; >= for 
every I — a/b. If r < 0, the system is in the in- 
sulating state < ipi >^ for at least one I. In the 
insulating or supersolid states, there must exist some 
kinds of charge density wave (CDW) or valence bond 
solid (VBS) orders which may be stabilized by longer 
range interactions or possible ring exch ang e interactions 
in Eqn[T] Up to an unknown prefactor [la ] , we can iden- 
tify the boson ( or vacancy ) densities on sites A and 
B and the boson kinetic energy on the link between A 
and B which are the order parameters for the CDW 
and VBS respectively as pa = i/j^ipa, Pb = i^i^b an d 
K AB = e^ s ^ b + e-^-^Va where Q = 2tt/3(1, 1) 
and x stands for dual lattice points only. We assume 
r < in Eqn[2j so the system is in the insulating state. 
In the following, we discuss the Ising limit first, then the 
easy-plane limit. 

4. Phase diagram (a) Ising limit. If 71 > 0, the sys- 
tem is in the Ising limit, the mean field solution is ip a = 
1, ipb = 0. The system is in the CDW order which could 
take checkboard (tt,tv) order [lij ]. Eqnl2]is an expansion 
around the uniform saddle point < V x A >= f = 1/2 
which holds in the SF and the VBS ( to be discussed 
in section 5 ). In the CDW state, a different saddle 
point where < VxA' >= 1 — a for sublattice A and 
< V x A b >~ a for sublattice B should be used. So the 
transition from the SF to the CDW is a strong first order 
transition. It can be shown [la] that there is only one 
vortex minimum ipb in such a staggered dual magnetic 
field with a < 1/2, the effective action inside the CDW 
state is: 



C 



CDW 



\(d, l -iA b M 2 +r\iJj b \ 2 +u\M 4 + 



+ ^(e^xd„A\-2ir5f6^) 2 



(3) 



where the vortices in the phase win ding of ij) b should be 
interpreted as the the boson number [17| . The gauge field 
A a is always massive. 

EqnELhas the structure identical to the conventional 
imponent Ginzburg-Landau model for a " super- 
conductor " in a "magnetic" field. By a duality trans- 
forma1(i2>)i back to the boson ^ picture, it can be shown 
that the transition driven by Sf is in the same universal- 
ity class of Mott to superfluid transition which has the 
exact exponents z = 2, v = 1/2 J? = with logarith- 
mic corrections first discussed in (lO, [Till- This fact was 
used in [18(] to show that for type II superconductors, 
the gauge field fluctuations will render the vortex fluid 
phase intruding at H c \ between the Messiner and the 
mixed phase Fig2a. For parameters appropriate to the 
cuprate superconductors, this intrusion occurs over too 
narrow an interval of H to be observed in experiments. 
In the present boson problem with the nearest neighbor 
interaction V\ > in Eqn[T] which stabilizes the (it, ir) 
CDW state at / = 1/2, this corresponds to a CDW su- 
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FIG. 1: (a) Bosons at filling factor / are hopping on a honey- 
comb lattice ( dashed line ) which has two sublattices A and 
B. Its dual lattice is a triangular lattice ( solid line ) which has 
three sublattices X,Y, Z. In the easy-plane limit, one of the 
three VBS states [14| is shown in the figure. The thick dashed 
bond is twice as strong as that of the thin dashed bond [l5| . 
The other two VBS can be obtained by R^tt/s or (b) 
and (c) are the zero temperature phase diagrams of the chemi- 
cal potential p versus r in Eqn[2]in the honeycomb lattice, (b) 
The Ising limit 71 > 0. There is a CDW supersolid (CDW-SS) 
state intruding between the commensurate CDW ( C-CDW ) 
state at / = 1/2 and the in-commensurate CDW (IC-CDW) 
state at 1/2 + 8f. The CDW-SS has the same lattice symme- 
try breaking as the C-CDW . (c) The Easy-Plane limit 71 < 0. 
There is a Valence Bond Supersolid (VB-SS) state intruding 
between the commensurate VBS (C-VBS) state at / = 1/2 
and the in-commensurate VBS ( IC-VBS ) state at 1/2 + Sf. 
The VB-SS has the same lattice symmetry breaking as the C- 
VBS . The thin ( thick ) line is the 2nd ( 1st ) order transition. 
The 1st order transition in the Ising ( Easy-plane ) limit is 
weakly ( strongly ) one. 



persolid (CDW-SS) state intruding between the commen- 
surate CDW state at / = 1/2 and the in-commensurate 
CDW state at 1/2 + 5/ which could be stabilized by fur- 
ther neighbor interactions in EqnQ] as shown in Fig. 2a. 
The hrst transition is in the z = 2, v = 1/2, r] = univer- 
sality class, while the second is a 1st order transition. We 
expect the intruding window at our q — 2 system is much 
wider than that of the q — 1 system. In the CDW-SS 
state, < ipb >= 0, there is the gapless superfluid mode 
represented by the dual gauge field A\ in EqnJ21 there is 
also the same (it, ir) diagonal order as the C-CDW where 
both A a and A^ are massive. For example, near the C- 
CDW to the CDW-SS transition, the superfluid density 



should scale as p s 



1/2| 



(d+z-2)v _ 



1/2 1 = Sf 



with logarithmic corrections. It is known the SF is sta- 
ble against changing the chemical potential ( or adding 
bosons ) in Fig. lb. There must be a transition from the 
CDW-SS to the SF inside the window driven by the quan- 
tum fluctuation r in the Fig. lb. The universality class 
of this transition is likely to be first order and will be 
investigated in a future publication [l5j |. 

5. Phase diagram (b) Easy-plane limit. If 71 < 
0, the system is in the easy-plane limit. This limit 
could be reached by possible ring exchange interactions 
Q in EqnQ] At q = 2, the mean field solution is 
t/; a = e l6a ,ipb = e l9b . pa = Pb = 1, so the two sub- 
lattices remain equivalent. The uniform saddle point 



< V x A >= / = 1/2 holds in both the SF and the 
VBS, so the transition from the SF to the VBS is a weak 
first order transition. The system has a VBS order, the 
kinetic energy K ab — cos(Q-x + 6L) where 6L = 6 a — 6b- 
Upto the quartic order, the relative phase 6L is undeter- 
mined. Higher order terms are needed to determine the 
relative phase. It was shown in [l5| that there are only 3 
sixth order invariants: C\ = \tp a \ 6 + \?Pb\ 6 , C2 = (IV^I 2 + 

|^| 2 )IV>a| 2 W 2 ,C 3 = (V:V>b) 3 + W^a) 3 = ACOS30_. 

Obviously, only the last term C3 can fix the relative 
phase. It is easy to show that both signs of A are equiva- 
lent, in sharp contrast to the square lattice where the two 
signs of C sq — A cos A6 lead to either Columnar dimer or 
plaquette pattern. One VBS with A > 0, 6L — tt was 
shown in Fig. la. The other two VBS states can be ob- 
tained by P^l%- I n contrast to the square lattice where 
C sq is irrelevant near the QCP, C3 may be relevant, so 
the transition from the SF to the VBS in Fig. lb could be 
a 1st order transition. Slightly away from the half-filling, 
Eqnj2] becomes: 



C 



VBS 



1 

4e 2 

2A cos 36» 



where 6± = a ±0b- 

Obviously, the sector is massive ( namely, 9 a and 9b 
are locked together ) and can be integrated out. Assum- 
ing A > 0, then 6*_ = it. Setting ip + = e l9 + ~ ip a ~ — ^ 
in Eqnj2] leads to Eqnj3] with u = 270, so the discus- 
sions on Ising limit case following Eqn[3] also apply. In 
the present boson problem with possible ring exchange 
interactions in EqnQ] which stabilizes the VBS state at 
/ = 1/2, this corresponds to a VBS supersolid (VB- 
SS) state intruding between the commensurate VBS ( 
C-VBS) state at / = 1/2 and the in-commensurate VBS 
(IC-VBS) state at 1/2+5/ as shown in Fig.2b. In this IC- 
VBS state, Sf valence bonds shown in Fig. la is slightly 
stronger than the others. In the VB-SS state, < ip a >=< 

1p b >= 0, but < Vl^a > = < Ipllpb >=- < 0, 

so there is a VBS order Kab — cos(Q • x + 6L) which is 
the same as the C-VBS, the superfluid density p s ~ 8f. 
Again, the first transition is in the z = 2,i/=l/2,?7 = 
universality class, while the second is 1st order. The na- 
ture of the transition from the VB-SS to the SF where 
< ipa >=< ipb >= and < ipl^pb >= inside the win- 
dow driven by the quantum fluctuation r in the Fig. 2b 
will be studied in [la |. 

6. Square lattice. As said below EqnJ^l with cor- 
respondingly defined ip a /b m a square lattice, upto the 
quartic level, Eqnl2]is the same as that in the square lat- 
tice derived in [1 31 ] - So in the Ising limit, Fig. 2b remains 
the same as Fig. lb. However, in the Easy-plane limit, 
as shown in [l3j , the lowest order term coupling the two 
phases 9 a / b is C sq — A cos AO. If A is positive ( negative), 
the VBS is Columnar dimer ( plaquette ) pattern. So the 
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FIG. 2: (a) The phase diagram slightly away from q — 1 
is the same as type-II superconductor in external magnetic 
field H where there should be a vortex liquid state intrud- 
ing between the Messiner state and the vortex lattice state. 
But the intruding regime is too narrow to be seen in type-II 
superconductors 18). (b) and (c) are the zero temperature 
phase diagrams of the chemical potential u versus r in Eqn[3] 
in square lattice, (b) The Ising limit 71 > is similar to 
Fig. lb. (c) The Easy-Plane limit 71 < is similar to Fig.lc 
except the SF to the C-VBS and the SF to the VB-SS could 
be 2nd order due to the so called "deconfined quantum critical 
point" . (b) and (c) are drawn only above half filling, u — > —fi 
corresponds to below half filling. The IC-CDW in (b) can be 
stabilized only by very long range interactions in Eqn[T] If it 
is not stable, then Fig.2b reduces to Fig. 3b. The first order 
transition in (b) is a strong one. 



C3 term in Eqn|4] need to be replaced by C sq . Because 
C sq is irrelevant near the QCP, the transition from the 
SF to the VBS could be a 2nd order transition through 
the deconfined QCP as shown in Fig. 2c. We expect 
that the SF to the VB-SS transition could also be 2nd 
order through a novel deconfined quantum critical line 
shown in Fig. 2c and will be studied in [lij]. 

7. Implication on Quantum Monte-Carlo (QMC) sim- 
ulations. QMC simulations of hard core bosons on square 
lattice with V\ and V% interactions find a stable striped 
0,0) and (0,tt) SS ( Fig.3b ) @. A stable (7r,7r) SS 
[8| can be realized in soft core boson case. But the 
nature of the CDW to supersolid transition has never 
been addressed. Our results show that the CDW to 
the SS transition must be in the same universality class 
of Mott to superfluid transition with exact exponents 
z = 2, v = l/2,r/ = with logarithmic corrections. 
It is important to (1) confirm this prediction by finite 
size scaling through the QMC simulations in square lat- 
tice for (0, ir) and (n, 0) supersolid in hard core case 
(Fig. 3b) and (7r,7r) supersolid in the soft core case (2) 
do similar things in honeycomb lattice to confirm Fig.l. 
(3) To Eqb[T] with U = 00, Vi > 0, adding ring ex- 
change term -K s Y^ijkli^i^k^ + h - c -) where i, j, k, I 
label 4 corners of a square in the square lattice and 
~ K h Y.i ki mn ( h \ h 3 h k b i h ii h n + h-c.) where k, I, m, n la- 
bel 6 corners of a hexagon in the honeycomb lattice to 
stabilize the C-VBS state at half filling ( K Sl Kh > are 
free of sign probelm in QMC ), then confirm the predic- 
tion on C-VBS to VB-SS transition in Fig. 2c and Fig.lc. 
The second transition ( CDW-SS to IC-CDW in Fig. lb 
or 2b and the VB-SS to IC-VBS in Fig lc or Fig.2c ) is 
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FIG. 3: ( The thin (thick) line is a 2nd (1st ) order transition, 
(a) Phase diagram for V% — 0, U — 00, the (7r,7r) SS, having 
negative compressibility, is unstable against phase separation 
into the SF and solid 7|. The DVM does not apply in this 
case, (b) Phase diagram for Vi = 0, U — 00, there is a narrow 
window of (tt, 0) SS sandwiched between the (ir, 0) solid and 
the SF. The SS disappears at half-filling. The DVM does 
apply in this case. It corresponds to Fig. 2b from the DVM. 



hard to be tested in QMC, because some very long range 
interactions are needed to stabilize the IC-CDW or the 
IC-VBS state. They are first order transition anyway. 

In fact, one of the predictions in this letter on the scal- 
ing of the superfluid density p s ~ \p — 1/2| was already 
found in the striped (7r, 0) solid to striped supersolid tran- 
sition by QMC in Sec.V-B [7]. In fact, as shown in section 
4, there should be logarithmic correction to the scaling 
of p s , it remains a challenge to detect the logarithmic 
correction in QMC. Of course, the superfluid density is 
anisotropic p x s > pi in the (n, 0) solid, but they scale in 
the same way with different coefficients [7( . Although the 
authors in [7] suggested it is a 2nd order transition, they 
did not address the universality class of the transition. 

8. Summary We studied some phases and phase tran- 
sitions in an extended boson Hubbard model near half 
filling on bipartite lattices such as honeycomb and square 
lattice. We identified boson density and boson kinetic 
energy operators to characterize symmetry breaking pat- 
terns in the insulating states and supersolid states. We 
found that in the insulating side, the transition at zero 
temperature driven by the chemical potential must be a 
C-CDW ( or C-VBS ) at half filling to a narrow window 
of CBW- ( VB-) supersolid, then to a IC-CDW ( IC-VBS 
) transition in the Ising ( easy-plane ) limit. The valence 
bond supersolid is a new kind of supersolid first proposed 
in this letter. The first transition is in the same univer- 
sality class as that from a Mott insulator to a superfluid 
driven by a chemical potential, therefore have exact expo- 
nents z = 2, v = 1/2, 77 = with logarithmic corrections. 
The second is a 1st order transition. The results achieved 
in this letter could guide QMC simulations to search for 
all these phases and confirm the universality class of the 
transitions. These transitions could be easily realized in 
near future atomic experiments in optical lattices. 
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